AD-A211  395 


fHF  (Wv 

CHEMICAL 
RESEARCH. 
DEVELOPMENT  & 
ENGINEERING 
CENTER 


CRDEC-CR-030 


IMPROVED  PATTERN  RECOGNITION 
TECHNIQUES  FOR  REAL  TIME  OPERATION 


Gary  W.  Small 

UNIVERSITY  OF  IOWA 
Iowa  City.  IA  52242 


May  1989 


DT1C 

ELECTE 
AUG  14 1989 


r  "f' Cf  rior-axiinm*  fcuK  1 

*  mkH  IllniM  Mi  ■ 
to  «albait«4» 


U-f.AH MY 
ARMAMENT 
MUNTTIONS 
CHEMICAL  COMMAND 


Abardaan  Proving  Ground.  Maryland  21010*5423 


8  ■  14  055 


Di  sclaimer 

The  findings  in  this  report  are  not  to  be  construed  as  an  official 
Department  of  the  Army  position  unless  so  designated  by  other 
authorizing  documents. 


Distribution  Statement 

Approved  for  public  release;  distribution  is  unlimited. 


UNCLASSIFIED 


URITY  CLASSIFICATION  OF  THIS  PAGE 


REPORT  DOCUMENTATION  PAGE 


Form  Approved 
OMB  No  0704-01 8B 


la.  REPORT  SECURITY  CLASSIFICATION 

UNCLASSIFIED 


2a.  SECURITY  CLASSIFICATION  AUTHORITY 


2b.  DECLASSIFICATION /DOWNGRADING  SCHEDULE 


4.  PERFORMING  ORGANIZATION  REPORT  NUMBER(S) 

CRDEC-CR-030 


1b  RESTRICTIVE  MARKINGS 


3  DISTRIBUTION /AVAILABILITY  OF  REPORT 

Approved  for  public  release; 
distribution  is  unlimited. 


5.  MONITORING  ORGANIZATION  REPORT  NUMBER(S) 


6a.  NAME  OF  PERFORMING  ORGANIZATION  6b.  OFFICE  SYMBOL  7a.  NAME  OF  MONITORING  ORGANIZATION 

(If  applicable) 

University  of  Iowa 


6c.  ADDRESS  (City,  State,  and  ZIP  Code) 


Iowa  City,  I A  52242 


8a.  NAME  OF  FUNDING  /SPONSORING 
ORGANIZATION 

CRDEC 


8c  ADDRESS  (City,  State,  and  ZIP  Code) 


8b.  OFFICE  SYMBOL 
(If  applicable) 
SMCCR-RSL-S 


7b.  ADDRESS  (C/fy,  State,  and  ZIP  Code) 


9  PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUMBER 

DAAA15-86-C-0034 


10.  SOURCE  OF  FUNDING  NUMBERS 


PROGRAM 
ELEMENT  NO. 


Aberdeen  Proving  Ground,  MD  21010-5423 


1 1 .  TITLE  (Include  Security  Classification) 

Improved  Pattern  Recognition  Techniques  for  Real  Time  Operation 


PROJECT 

TASK 

NO 

NO. 

12.  PERSONAL  AUTHORS) 

Smal 1 ,  Gary  W. 


13a.  TYPE  OF  REPORT 

Contractor 


16.  SUPPLEMENTARY  NOTATION 


13b  TIME  COVERED  114.  DATE  OF  REPORT  (Year,  Month.  Day)  115.  PAGE  COUNT 

from  86  Apr  to  88  Aug  I  1989  May  I  49 


GROUP  SUB-GROUP 


COR:  Robert  Kroutil,  SMCCR-RSL-S,  (301)  671-3021 


17.  COSATI  CODES  |  ■L^SUBJECT  TERMS  (Continue  on  reverse  if  necessary  and  identify  by  block  number) 

^^ourier  spectroscopy  >  \ 

\  ^Automatic  spectral  interpretation  __  ^ 

| .  .VStandoff  chemical  vapor  detection  al  gorithm.  i'f- 

BSTRACT  (Continue  on  reverse  if  necessary  and  identify  by  block  number)  _ 

Passive  Fourier  transform  infrared  remote  sensors  have  shown  tremendous  potential  for  use 
in  the  automated  detection  of  atmospheric  pollutants.  To  be  maximally  effective  as  an 
early  warning  alarm  device,  however,  passive  infrared  sensors  must  be  able  to  operate  in 
a  variety  of  infrared  backgrounds  and  under  conditions  in  which  the  background  is  rapidly 
changing.  For  economic  reasons,  the  sensor  should  be  the  simplest  device  that  can  meet 
the  aforementioned  operational  requirements.  In  the  work  described  here,  new  signal 
processing  strategies  have  been  developed  that  allow  target  chemical  species  to  be  detec¬ 
ted  directly  from  time-domain  infrared  data.  This  methodology  eliminates  the  Fourier 
transform  computation  through  the  development  of  digital  filters  keyed  to  the  character¬ 
istic  infrared  frequencies  of  the  target  molecules.  Through  this  approach,  both  the 
computational  and  instrumental  requirements  for  the  sensor  are  significantly  reduced. 

The  performance  of  this  methodology  will  be  illustrated  with  a  variety  of  test  data. 


VV.J  l  |\U44  ,1  T 


■-j  “5u\ H\V 


20  DISTRIBUTION /AVAILABILITY  OF  ABSTRACT 
E  uNCLASSltlED/UNUMITED  □  SAME  AS  RPT 
22a.  NAME  OF  RESPONSIBLE  INDIVIDUAL 

SANDRA  J.  JOHNSON 


DO  Form  1473,  JUN  86  / 


21  ABSTRACT  SECURITY  CLASSIFICATION 

□  dtic  users  UNCLASSIFIED 


22b  TELEPHONE  (Include  Area  Code)  22c  OFFICE  SYMBOL 

(3011  671-2914  SMCCR-SPS-T 


Previous  editions  are  obsolete. 


SECURITY  CLASSIFICATION  OF  THIS  PAGE 


UNCLASSIFIED 


UNCLASSIFIED 


2 


UNCLASSIFIED 


PREFACE 


The  work  described  in  this  report  was  authorized  under  Contract  No. 
DAAA15-86-C-0034.  This  work  was  started  in  April  1986  and  completed  in 
August  1988. 

The  use  of  trade  names  or  manufacturers'  names  in  this  report  does 
not  constitute  an  official  endorsement  of  any  commercial  products.  This 
report  may  not  be  cited  for  purposes  of  advertisement. 

Reproduction  of  this  document  in  whole  or  in  part  is  prohibited 
except  with  permission  of  the  Commander,  U.S.  Army  Chemical  Research, 

Development  and  Engineering  Center,  ATTN:  SMCCR-SPS-T,  Aberdeen  Proving  Ground, 
Maryland  21010-5423.  However,  the  Defense  Technical  Information  Center  and  the 

National  Technical  Information  Service  are  authorized  to  reproduce  the  document 
for  U.S.  Government  purposes. 

This  document  has  been  approved  for  release  to  the  public. 

ACKNOWLEDGMENTS 

The  work  described  in  this  report  represents  the  combined  efforts  of 
a  group  of  analytical  chemists  working  in  the  Department  of  Chemistry  at  the 
University  of  Iowa.  The  contributions  made  by  Abigail  Barber,  Joanne  Bjerga, 
Scott  Carpenter,  Amy  Harms,  and  Thomas  Kaltenbach  to  the  work  are  gratefully 
acknowledged.  Dr.  Robert  T.  Kroutil  and  coworkers  of  CRDEC  are  also  acknow¬ 
ledged  for  supplying  the  passive  FTIR  data  used  in  the  work  and  for  their 
participation  in  many  invaluable  technical  discussions. 


3 


Blank 


4 


CONTENTS 


INTRODUCTION  .  9 

EXPERIMENTAL  SECTION  .  9 

RESULTS  AND  DISCUSSION  .  10 

Overview  of  Passive  FTIR  Methodology  and  Conventional 

Data  Processing  . 10 

Motivation  for  Alternative  Algorithm  Development  Work  .  12 

Development  of  Digital  Filtering  Methodology  .  18 

Mathematical  Background  .  18 

Regression  Approach  to  FIR  Filter  Generation  .  19 

Generation  of  FIR  Filters  by  Stepwise  Regression  .  20 

FIR  Matrix  Filters .  23 

Quantitative  Comparison  of  SF§  Filter  Performance  .  32 

Quantitative  Nature  of  Filter  Responses  . 34 

Pattern  Recognition  Methodology  Based  on  Filtered 

Interferogram  Segments . . .  35 

Overview  of  Methodology  . 35 

Construction  of  Linear  Discriminants . 38 

Testing  of  Computed  Discriminant  .  39 

Critical  Evaluation  of  Time-Domain  Methodology  .  43 

Questions  of  Spectral  Resolution  . 43 

Effects  of  Phase  Errors  .  43 

Source  of  False  Alarms . 44 

CONCLUSIONS  .  46 

RECOMMENDATIONS  .  47 

REFERENCES .  49 


5 


LIST  OF  FIGURES 


Page 


Figure 
No. 

1.  A  1024-point  Interferogram  collected  by  a  passive  FTIR' 

spectrometer  mounted  on  a  helicopter . .  11 

2.  The  spectrum  arising  from  the  application  of  the  FFT  to  the 

Interferogram  depicted  in  Figure  1  .  11 

3.  Difference  spectrum  resulting  from  the  subtraction  of  a  background 

reference  spectrum  from  the  sample  spectrum  of  Figure  2  .  13 

4.  A  typical  laboratory  difference  spectrum  showing  the  absorption 

bands  of  butyric  acid  superimposed  on  a  flat  baseline  .  13 

5.  Segments  from  six  passive  FTIR  Interferograms  .  15 

6.  The  region  from  900-1000  cm"!  Is  plotted  In  the  transformed  spectra 
corresponding  to  the  Interferograms  from  which  the  segments  In 

Figure  5  were  drawn  .  15 

7.  The  frequency-dependent  action  of  a  digital  filter  . .  16 

8.  Three  spectra  are  plotted  from  650-1250  cnfl  . 16 

9.  Points  1-50  In  the  Interferograms  corresponding  to  the  three 

spectra  In  Figure  8  .  17 

10.  Points  50-150  In  the  Interferograms  corresponding  to  the  three 

spectra  In  Figure  8  .  17 

11.  Points  150-250  In  the  Interferograms  corresponding  to  the  three 

spectra  In  Figure  8  . . 18 

12.  Points  175-250  In  the  six  Interferogram  segments  of  Figure  5  after 

applying  the  20-term  filter  computed  by  stepwise  regression . .  21 

13.  Sum  of  squares  after  applying  the  20-term  filter  vs.  Interferogram 

number  for  each  Interferogram  In  a  typical  helicopter  data  run  .  21 

14.  Sum  of  squares  after  applying  the  20-term  filter  vs.  Interferogram 
number  for  each  Interferogram  In  a  typical  ground  vehicle  data  run  .  22 

15.  Sum  of  squares  vs.  Interferogram  number  for  the  ideal  response  of 

a  Gaussian  filter  centered  at  941.1  cm”l  and  having  a  FWHM  of  54.0 
cm~^  . 23 

16.  Plot  of  percentage  of  variance  explained  vs.  Interferogram  point  ...  25 

17.  Plot  of  number  of  filter  coefficients  required  vs.  Interferogram 

point  .  25 

18.  Plot  of  percentage  variance  explained  In  the  matrix  filter 

generation  vs.  intensity  of  filtered  Interferogram  points  .  26 


6 


19. 


Plot  of  points  175-250  in  six  interferogram  segments  after 
application  of  the  matrix  filter  . 


26 


20.  Sum  of  squares  over  points  175-250  after  applying  the  matrix 

filter  vs.  Interferogram  number  for  the  helicopter-based  data  run  ..  28 


21.  Sum  of  squares  over  points  175-250  after  applying  the  matrix 

filter  vs.  Interferogram  number  for  the  ground-based  data  run  .  28 

22.  Sum  of  squares  over  51  points  exceeding  the  9055  variance  cutoff 

vs.  Interferogram  number  for  the  helicopter-based  data  run  .  29 

23.  Sum  of  squares  over  51  points  exceeding  the  9055  variance  cutoff 

vs.  Interferogram  number  for  the  ground-based  data  run  .  29 

24.  Frequency  response  function  for  the  DMMP  matrix  fll+er  .  30 

25.  Sum  of  squares  over  points  175-250  vs.  Interferogram  number  for 

the  Ideal  response  of  the  DMMP  filter . . . . . 30 

26.  Plot  of  percentage  of  variance  In  the  dependent  variable  explained 

vs.  Interferogram  point  for  the  DMMP  matrix  filter  .  31 

27.  Plot  of  number  of  filter  terms  vs.  Interferogram  point  for  the 

DMMP  matrix  filter . . . . . 31 

28.  Sum  of  squares  over  47  points  exceeding  the  8055  variance  cutoff 

vs.  Interferogram  number  for  the  DMMP  data  run  . 32 

29.  Points  175-250  In  six  Interferograms  derived  from  the  helicopter- 

based  data  run  at  1000  ft  . . .  36 

30.  The  region  from  920-970  cm-*  in  the  spectra  derived  from  the 

Interferograms  In  Figure  29  .  36 

31.  Points  175-250  In  the  six  Interferograms  after  application  of  the 

SFg  matrix  filter . 37 

32.  Plot  of  the  magnitude  of  the  filtered  Interferogram  segment  vs. 

Integrated  area  of  the  SFg  spectral  band  . . . .  37 

33.  Two-point  mean  discriminant  score  vs.  interferogram  number  for  a 

ground-based  data  run  . .  40 

34.  Two-point  mean  discriminant  score  vs.  Interferogram  number  for  a 

helicopter-based  data  run  . .  40 

35.  Two-point  mean  discriminant  score  vs.  Interferogram  number  for  a 

helicopter-based  data  run  . . 41 

36.  Two-point  mean  discriminant  score  vs.  Interferogram  number  for  a 

ground-based  data  run  .  41 

37.  Plot  of  percentage  variance  explained  in  the  matrix  filter 
calculation  vs.  computed  discriminant  value  for  Interferogram 

points  175-250  . 42 


7 


_>d.  Transformed  spectra  corresponding  to  a  76-point  segment  taken 
from  an  Interferogram  composed  of  a  single  cosine  frequency 
(941.1  cm"1)  .  44 

39.  Two-point  mean  discriminant  score  plot  for  the  helicopter-based 

data  run  after  Introduction  of  catastrophic  phase  errors  .  45 

40.  Two-point  mean  discriminant  score  plot  for  the  ground-based  data 

run  after  Introduction  of  catastrophic  phase  errors  .  45 

41.  Transformed  spectrum  of  Interferogram  503  In  Figure  35  . .  46 

LIST  OF  TABLES 

Table  No.  Page 

I  Slgnal-to-Nolse  Comparison  of  Filter  Performance  SFg 

Hel  1  copter-Based  Oata  .  33 

II  Slgnal-to-Nolse  Comparison  of  Filter  Performance  SFg 

Ground-Based  Data . . .  33 

III  Comparison  of  Computational  Efficiency  . .  34 


8 


IMPROVED  PATTERN  RECOGNITION  TECHNIQUES  FOR  REAL  TIME  OPERATION 


INTRODUCTION 


Passive  Infrared  spectrometers  have  come  Into  recent  use  for  remote  sensing 
of  atmospheric  pollutants!-!®.  These  devices  consist  of  a  Fourier  transform 
infrared  (FTIR)  spectrometer  In  which  the  Infrared  source  has  been  removed  and 
a  set  of  specialized  Input  optics  added.  The  use  of  passive  FTIR  Instruments 
Is  based  on  the  premise  that  most  gaseous  species  exhibit  a  spectral  signature 
In  the  Infrared  region  that  will  be  superimposed  on  the  normal  background 
Infrared  emission  In  the  fleld-o f-view  of  the  spectrometer.  When  controlled 
by  a  dedicated  computer  and  appropriate  software#  the  Instrument  can  serve  as 
an  automated  alarm  system  for  the  detection  of  specific  target  gases. 

Applications  for  passive  FTIR  Instruments  have  Included  both  stationary 
and  mobile  environments.  Recently#  Increased  attention  has  been  given  to  the 
mobile  application  for  the  following  reasons:  (1)  the  Instrument  Is  allowed 
to  operate  In  the  cleaner  environment  of  a  ground  or  airborne  vehicle;  (2)  a 
single  Instrument  can  be  used  to  monitor  a  larger  area;  and  (3)  with  fewer 
Instruments  needed#  more  highly  trained  personnel  can  be  used  to  oversee  the 
Installation#  maintenance#  and  operation  of  the  system. 

In  a  mobile  application#  the  spectrometer  must  be  able  to  operate  against 
a  variety  of  Infrared  spectral  backgrounds  and  must  be  able  to  resist  the 
effects  of  rapid  background  changes.  In  a  high-speed  vehicle  (e.g.  helicopter 
or  other  aircraft)#  the  data  processing  tasks  must  be  accomplished  rapidly. 

This  Includes  the  complete  decision-making  procedure  for  the  determination  of 
the  presence  or  absence  of  the  target  analyte(s).  Additionally#  for  the 
system  to  be  practical,  the  infrared  Instrumentation  must  be  both  reliable  and 
economical.  Both  reliability  and  cost  are  Inversely  related  to  the  degree  of 
sophistication  required  In  both  the  Instrument  and  the  controlling  computer. 
Further#  the  required  sophistication  of  both  the  Infrared  Instrumentation  and 
the  computational  hardware  Is  directly  related  to  the  amount  of  data  that  must 
be  collected  and  the  requirements  for  processing  that  data. 

This  report  describes  efforts  at  the  University  of  Iowa  to  reduce  both 
Instrumental  and  computational  requirements  for  passive  Infrared  sensors 
through  the  development  of  novel  data  processing  algorithms.  In  addition# 
algorithms  are  being  sought  that  are  resistant  to  changes  In  the  Infrared 
background.  These  algorithms  are  described  In  detail  and  evaluated  with  a 
variety  of  test  data. 


EXPERIMENTAL  SECTION 


A  passive  Infrared  Interferometer  was  constructed  by  Honeywell 
Corporation  to  specifications  provided  by  the  U.S.  Army  Chemical  Research# 
Development#  and  Engineering  Center#  Edgewood#  MD.  The  system  consisted  of  a 
flex-pivot  "porch  swing"  Mlchelson  Interferometer  based  on  the  design  of 
Walker  and  Rex!!,  j^e  spectrometer  used  a  liquid-nitrogen  cooled  Hg:Cd:Te 
detector  to  collect  spectral  background  radiation  from  the  8  to  12  micron 
region.  For  the  work  reported  here#  1024-point  Interferograms  were  collected# 
with  a  corresponding  spectral  resolution  of  approximately  4  cm-!,  The 
standard  data  collection  rate  for  this  instrument  was  5  Interferograms/ sec. 

In  order  to  test  the  feasibility  of  a  mobile  application  for  the  sensor 
and  to  collect  data  under  a  variety  of  background  conditions#  the  spectrometer 
was  mounted  on  a  shock-absorbing  platform  and  placed  In  both  a  helicopter  and 
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a  ground  vehicle.  Several  runs  of  data  were  collected  at  a  variety  of  speeds 
and  distances  from  ground  sources  of  SFg  and  dimethyl  methyl phospnonate 
(DMMP).  These  gases  were  used  to  simulate  the  presence  of  target  analytes 
that  should  be  detected  by  the  spectrometer.  SFg  has  a  single  strong 
absorption  at  940  cm“l,  while  DMMP  has  characteristic  bands  at  820,  910,  and 
1040  cm"l.  The  data  collection  was  performed  at  Edgewood,  MD  and  Dugway,  UT. 
The  collected  Interferograms  were  written  onto  magnetic  tape  and  shipped  to 
the  University  of  Iowa  where  the  algorithm  development  work  was  performed. 

All  computer  software  used  In  this  work  was  written  In  FORTRAN  77  and 
Implemented  on  a  PRIME  9955  Interactive  computer  system  operating  at  the 
Gerard  P.  Weeg  Computing  Center  at  the  University  of  Iowa.  Some  computations 
made  use  of  subroutines  from  the  IMSL  statistical  package  (IMSL,  Inc., 
Houston,  TX).  Some  graphics  output  made  use  of  the  TELAGRAF  Interactive 
graphics  package  (Integrated  Software  Systems  Corp.,  San  Diego,  CA) .  A 
Hewlett-Packard  7475A  digital  plotter  was  used  as  the  output  device. 


RESULTS  AND  DISCUSSION 


flx.eo.lea  el  Passive  ftir  Matti.Qdola.gy  end  flon volitional  Ante  P.roc.esslng 

In  a  passive  FTIR  spectrometer,  the  initial  Infrared  data  Is  collected  In 
the  form  of  an  1 nterferogram,  a  composite  waveform  composed  of  periodic 
signals  corresponding  to  each  of  the  detected  Infrared  frequencies.  In 
conventional  spectral  processing,  a  Fast  Fourier  Transform  (FFT)  is  applied  to 
the  Interferogram  to  decompose  It,  thereby  obtaining  the  amplitudes  of  the 
Individual  frequency  components  present.  Taken  together,  these  amplitudes 
define  the  spectrum  of  Infrared  light  being  detected  by  the  spectrometer.  The 
characteristic  Infrared  frequencies  of  any  chemical  compounds  present  are 
superimposed  on  this  overall  spectrum  as  either  absorption  or  emission  bands. 
Figure  1  depicts  a  typical  1024-point  FTIR  Interferogram,  and  Figure  2 
displays  the  corresponding  spectrum  yielded  by  the  FFT.  The  arrow  In  Figure  2 
points  to  the  characteristic  band  of  SFg,  arising  here  as  an  absorption  (l.e. 
a  decrease  In  light  Intensity). 

In  conventional  processing  algorithms,  the  detection  of  a  target  gas 
begins  with  the  collection  of  a  representative  background  emission 
Interferogram  and  an  Interferogram  collected  with  the  potential  toxic  gas 
source  In  the  fie Id-of-vlew  of  the  Instrument.  The  FFT  Is  then  applied  to  the 
collected  Interferograms.  The  difference  spectrum  resulting  from  subtracting 
the  sample  and  background  spectra  Is  analyzed  for  the  spectral  signature  of 
the  compound  being  monitored.  Under  conditions  In  which  the  background  Is 
changing,  however,  the  difference  spectrum  can  contain  anomalous  features, 
possibly  resulting  In  either  false  alarms  or  missed  alarms. 

As  an  Illustration,  the  Interferogram  depicted  In  Figure  1  was  collected 
when  the  spectrometer  was  mounted  In  a  helicopter  flying  at  25  knots  and  at  an 
altitude  of  approximately  200  feet.  A  reference  spectrum  was  obtained  by 
transforming  an  Interferogram  collected  two  seconds  earlier.  Based  on  the 
position  of  the  helicopter  and  an  Inspection  of  the  spectrum,  no  SFg  Is 
Indicated  In  this  reference  spectrum.  The  spectrum  in  Figure  2  (l.e.  the 
sample  spectrum  In  this  example)  and  the  reference  spectrum  were  each 
normalized  to  unit  area  to  correct  for  overall  differences  in  Infrared 
Intensity.  Figure  3  depicts  the  difference  spectrum  obtained  by  subtracting 
the  reference  spectrum  from  that  of  the  sample.  The  SFg  absorption  remains 
evident,  but  a  number  of  spectral  artifacts  have  been  Introduced  that  may 
Interfere  with  numerical  decision-making  algorithms  designed  to  produce  a  yes 
or  no  regarding  the  presence  of  SFg. 
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Figure  1.  A  1024-point  Interferogram  collected  by  a  passive  FTIR  spectrometer 
mounted  on  a  helicopter.  Light  Intensity  detected  Is  plotted  vs. 

Interferogram  point  number.  The  Interferogram  has  been  rotated  such  that  the 
maximum  Intensity  (center  burst)  Is  at  point  1.  A  ground  source  of  SFg  was  In 
the  fleld-of-vlew  of  the  Instrument  when  this  Interferogram  was  collected. 


gure  2.  The  spectrum  arising  from  the  application  of  the  FFT  to  the 
interferogram  depicted  In  Figure  1.  Light  Intensity  Is  plotted  vs.  spectral 
frequency  (cm-1).  The  absorption  band  due  to  the  presence  of  SFg  Is  Indicated 
by  the  arrow. 
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By  way  of  comparison,  the  results  of  a  typical  laboratory  spectral 
subtraction  are  depicted  In  Figure  4.  A  normalized  reference  spectrum  was 
subtracted  from  a  sample  spectrum  of  butyric  acid  to  produce  the  difference 
spectrum  shown  In  the  figure.  These  spectra  were  collected  on  an  IBM  M98 
laboratory  FTIR  spectrometer  under  stable  background  conditions.  The 
characteristic  absorption  bands  of  butyric  acid  are  clearly  evident  against  a 
flat  baseline.  This  example  clearly  demonstrates  the  severity  of  the 
artifacts  Introduced  In  Figure  3  by  an  unstable  spectral  background. 

A  variety  of  techniques  have  been  Investigated  to  help  overcome  the 
presence  of  these  artifacts.  Digital  filters  have  been  employed  to  attempt  to 
discriminate  against  the  artifacts.  In  addition,  attempts  have  been  made  to 
characterize  common  background  features  In  hopes  of  accounting  for  the 
presence  of  background  changes.  All  such  methods  are  Inherently  limited, 
however,  as  an  Infinite  variety  of  possible  background  changes  exist. 

Extensive  Inspections  of  spectra  collected  under  moving  background 
conditions  reveal  one  overriding  principles  spectral  subtraction  Is  a 
laboratory  procedure  based  on  the  assumption  that  a  stable  spectral  background 
can  be  established.  The  use  of  spectral  subtraction  techniques  under 
conditions  of  changing  backgrounds  Is  an  Invalid  application  of  the  technique. 

Motivation  ion  Alternative  Algorithm  Development  Work 

The  spectral  subtraction  procedure  described  above  Is  useful  In  a 
laboratory  environment  because  It  represents  a  conceptually  simple  way  to 
extract  an  absorption  or  emission  band  from  the  Infrared  background  response, 
provided  the  background  response  Is  known.  The  procedure  Is  applied  to 
spectra  because  the  FFT  represents  the  conceptually  simplest  way  to  remove 
Information  regarding  spectral  frequencies  not  pertinent  to  the  detection  of 
the  analyte. 

Often,  the  conceptually  simplest  approach  to  a  data  analysis  problem  Is 
not  the  most  computationally  efficient.  The  standard  Cooley-Tukey  FFT 
requires  over  10,000  multiplications  for  a  single  1024-point  Interferogram. 
Additionally,  as  argued  In  the  discussion  above,  this  conceptually  simplest 
approach  based  on  spectral  subtraction  Is  not  really  appropriate  for  the 
present  application.  The  goal  for  our  research  has  been  to  develop  a  data 
processing  algorithm  for  mobile  passive  FTIR  spectrometers  that  makes  better 
use  of  the  available  computational  resources,  while,  at  the  same  time, 
overcoming  the  problems  of  changing  spectral  backgrounds. 

As  motivated  above,  the  goal  for  such  an  algorithm  Is  twofold.  First, 
the  Information  pertaining  to  the  characteristic  frequencies  of  a  target 
molecule  must  be  separated  from  the  Information  pertaining  to  other 
frequencies.  Second,  any  absorption  or  emission  band  present  In  the  frequency 
window  of  Interest  must  be  extracted  from  the  spectral  baseline.  This 
extracted  Information  can  then  be  used  with  a  standard  decision-making 
(pattern  recognition)  algorithm  to  produce  a  yes/no  response  regarding  the 
presence  of  the  target  analyte. 

As  Information  regarding  the  characteristic  frequencies  of  a  target 
molecule  Is  also  present  In  the  raw  Interferogram,  it  can  be  argued  that  use 
of  the  Fourier  transform  Is  an  Inefficient  strategy  for  frequency  selection. 
The  problem  with  an  Interferogram-based  analysis,  however,  centers  on  the 
composite  nature  of  the  waveform.  The  periodic  signals  corresponding  to  the 
Individual  frequencies  co-add  to  form  the  Interferogram.  Information 
regarding  the  specific  frequencies  of  Interest  Is  therefore  obscured  by 
frequency  Information  that  has  no  bearing  on  determining  the  presence  of  the 
target  compound. 

This  problem  Is  Illustrated  by  an  Inspection  of  Figure  5.  Segments 
(points  175-250  In  the  1024-point  rotated  Interferogram)  from  six 
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Figure  3.  Difference  spectrum  resulting  from  the  subtraction  of  a  background 
reference  spectrum  from  the  sample  spectrum  of  Figure  2.  Difference  In 
Intensity  Is  plotted  vs.  frequency  (cm"1).  The  artifacts  Introduced  are  due 
to  changes  In  the  background  between  the  times  the  two  spectra  were  collected. 


Figure  4.  A  typical  laboratory  difference  spectrum  showing  the  absorption 
bands  of  butyric  acid  superimposed  on  a  flat  baseline.  Difference  In 
Intensity  Is  again  plotted  vs.  frequency  (cm"1).  This  type  of  spectrum 
results  when  a  stable  spectral  background  exists. 
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Interferograms  are  plotted.  Figure  6  depicts  the  region  of  900-1000  cm"*  in 
the  corresponding  transformed  and  normalized  spectra.  Counting  from  the  top 
of  the  figure#  spectra  3#  5#  and  6  appear  to  contain  SFg  information/  although 
the  bands  are  very  weak.  It  Is  apparent  from  an  inspection  of  Figure  5/ 
however/  that  it  would  be  difficult/  at  best/  to  base  a  detection  procedure  on 
the  raw  Interferogram  Information  alone. 

Time-domain  digital  f11ters*2  are  mathematical  transforms  that  operate  on 
time-domain  data  In  a  frequency-dependent  manner.  The  result  of  the 
transformation  Is  a  filtered  time-domain  signal  that  has  had  Information 
pertaining  to  certain  frequencies  suppressed.  The  frequency  dependence  of  a 
filter  is  encoded  In  its  frequency  response  function.  In  effect#  this 
function  Is  a  spectrum  that  defines  which  frequencies  the  filter  will  pass. 
Figure  7  demonstrates  the  action  of  such  a  filter  in  the  spectral  domain.  The 
upper  plot  In  the  figure  Is  the  same  passive  FTIR  spectrum  depicted  In  Figure 
2.  The  middle  plot  is  a  filter  frequency  response  function  centered  on  the 
SFg  absorption  band.  This  Is  a  Gausslan-shaped  function  with  a  full-width  at 
half-maximum  (FWHM)  of  71.1  cm"*.  The  result  of  applying  this  filter  to  the 
FTIR  spectrum  Is  depicted  at  the  bottom  of  the  figure.  The  filter  effectively 
extracts  the  SFg  band  from  the  spectral  baseline. 

The  above  discussion  focused  on  the  familiar  spectral  domain.  The 
advantage  of  a  time-domain  filter#  however#  Is  that  a  filter  with  the 
frequency  response  depicted  In  Figure  7  can  operate  directly  on  the 
Interferogram#  producing  a  filtered  Interferogram  whose  transformed  spectrum 
Is  Identical  to  that  shown  at  the  bottom  of  the  figure.  In  this  manner#  the 
Fourier  transform  has  been  avoided  and  the  spectral  band  of  Interest  has  been 
extracted  from  the  baseline. 

Figures  8-11  demonstrate  graphically  that  an  Interferogram-based 
detection  scheme  Is  a  feasible#  straightforward  approach  to  the  passive  FTIR 
detection  problem.  Figure  8  depicts#  from  top  to  bottom#  a  Gausslan-shaped 
filter  frequency  response#  the  spectrum  resulting  from  applying  this  filter  to 
a  synthetic  Lorentzlan  absorption  band#  and  the  spectrum  resulting  from 
applying  the  filter  to  a  synthetic  Lorentzlan  emission  band.  The  two 
Lorentzlan  bands  are  Identical  In  shape#  width#  and  area. 

If  the  Inverse  FFT  Is  applied  to  the  three  spectra#  the  corresponding 
Interferograms  are  obtained.  Figures  9-11  depict  points  1-50#  50-150,  and 
150-250#  respectively#  In  these  three  computed  Interferograms.  An  Inspection 
of  these  figures  reveals  three  fundamental  facts  that  allow  an  Interferogram- 
based  detection  scheme  to  work. 

First#  a  separation  of  Information  exists  In  the  Interferogram  based  on 
the  width  of  the  spectral  feature.  The  three  segments  in  Figure  9  are 
Identical.  This  Indicates  that  the  bulk  of  the  Information  pertaining  to  the 
Gaussian  filter  function  Is  contained  In  the  first  50  Interferogram  points. 

The  time-domain  form  of  the  Gaussian  has  effectively  damped  out  by  point  110. 
The  Information  from  points  110-250  pertains  almost  entirely  to  the  Lorentzlan 
band.  This  observation  Implies  that  only  a  short  segment  of  the  Interferogram 
may  be  necessary  In  order  to  detect  the  presence  of  a  particular  spectral 
band.  By  collecting  only  a  short  segment  (termed  windowing  the 
Interferogram)#  Instrumental  requirements  are  reduced  and  many  background 
artifacts  are  automatically  eliminated. 

Second#  the  only  difference  between  the  absorption  and  emission  spectra 
exists  In  the  region  of  points  50-100.  Beyond  point  100#  the  signals  are 
Identical  In  magnitude#  differing  only  In  phase.  Stated  differently#  the 
lower  two  plots  In  Figure  11  are  exact  mirror  Images.  This  observation 
confirms  that  absorption  and  emission  bands  can  be  detected  In  the 
Interferogram  with  equal  ease. 

Third#  Information  regarding  the  Intensity  of  the  absorption  or  emission 
band  Is  encoded  In  the  Intensity  of  the  periodic  Interferogram  signal  produced 
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Figure  5.  Segments  from  six  passive  FTIR  Interferograms.  Light  Intensity  Is 
plotted  vs.  Interferogram  point  for  points  175-250. 


Figure  6.  The  region  from  900-1000  cm-1  Is  plotted  In  the  transformed  spectra 
corresponding  to  the  Interferograms  from  which  the  segments  In  Figure  5  were 
drawn.  Small  SFg  bands  around  940  cm"1  are  evident  In  the  third,  fifth,  and 
sixth  spectra,  counting  from  the  top  of  the  figure. 
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Figure  7.  The  frequency-dependent  action  of  a  digital  filter.  From  top  to 
bottom:  a  typical  passive  FTIR  spectrum*  a  filter  frequency  response 
function,  the  result  of  applying  the  filter  to  the  spectrum. 


Figure  8.  Three  spectra  are  plotted  from  650-1250  cm"1.  Top:  A  Gaussian 
filter  frequency  response.  Middle:  The  result  of  applying  this  filter  to 
extract  a  Lorentzlan  absorption  band.  Bottom:  The  result  of  extracting  an 
Identical  Lorentzlan  emission  band. 
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Figure  11.  Points  150-250  in  the  Interferograms  corresponding  to  the  three 
spectra  In  Figure  8. 

by  the  band.  This  suggests  that  a  possible  qualitative  Indicator  of  the 
presence  of  an  absorption  or  emission  may  be  a  simple  power  calculation  (e.g. 
sum  of  squares  of  Intensity)  over  an  appropriate  Interferogram  segment. 
Moreover,  such  a  computed  value  may  have  quantitative  spectral  Information 
also. 

The  above  discussion  demonstrates  unequivocally  that  the  combination  of 
windowing  the  Interferogram  and  applying  a  time-domain  digital  filter  can  be 
used  to  detect  the  presence  of  specific  absorption  or  emission  bands.  Through 
this  procedure,  the  band  Information  can  be  extracted  from  the  Infrared 
background  without  the  need  for  the  Fourier  transform  or  for  any  type  of 
spectral  subtraction  or  background  characterization.  The  development  and 
application  of  this  methodology  to  passive  FTIR  data  will  now  be  discussed. 

Development  Sd  Digital  Filtering  Methodology 

Mathematical  Background.  In  the  discussion  above,  the  frequency-domain 
form  of  a  digital  filter  was  applied  to  a  spectrum  by  multiplying  the  spectrum 
by  the  frequency  response  function  of  the  filter.  This  calculation  has  a 
corresponding  form  In  the  time  domain  that  Is  expressed  by  the  convolution 
Integral 

00 

y(t)  *  /htk)x(t-k)  dk  =  H(f)X(f)  (1) 

-os 

where  y(t)  Is  the  filtered  Interferogram,  the  product,  H(f)X(f),  Is  the 
filtered  spectrum,  H(f)  Is  the  frequency  response  of  the  filter,  X(f)  Is  the 
transformed  spectrum,  x(t)  Is  the  raw  Interferogram,  and  h(t)  Is  the  Inverse 
Fourier  transform  of  the  frequency  response,  termed  the  Impulse  response  of 
the  filter.  The  terms,  H(f)  and  X(f),  are  functions  of  frequency,  f,  while 
the  terms,  y(t)»  x(t)»  and  h(t),  are  functions  of  the  time  variables,  t  or  k. 

The  above  equation  pertains  to  continuous  functions.  For  digitized  data, 
the  discrete  convolution  Integral  has  the  form 
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(2) 
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where  Is  one  point  on  the  filtered  Interferogram#  and  h^  and  x-^-k  define 
discrete  points  on  the  Impulse  response  function  and  Interferogram# 
respectively.  Given  that  the  interferogram  Is  only  sampled  over  a  finite 
range*  the  discrete  form  of  the  convolution  Integral  must  be  approximated  by 
use  of  a  finite  series  of  terms.  This  approximation  most  often  takes  the  form 

yt  =  h0xt  +  ...  +  hnx^_n  (3) 

where  n+1  terms  are  summed  to  estimate  each  filtered  interferogram  point#  y^. 
The  terms#  h|<#  can  be  considered  weighting  coefficients  that  determine  the 
frequency  dependence  of  the  filter.  Effectively#  each  filtered  Interferogram 
point  Is  formed  from  a  linear  combination  of  the  corresponding  raw  data  point 
and  a  set  of  preceding  raw  data  points.  This  type  of  filter#  defined  by  the 
h|<#  Is  termed  a  finite  Impulse  response  (FIR)  filter. 

The  hj<  terms  must  be  found  that  best  approximate  the  action  of  the  filter 
frequency  response  function.  This  Is  most  often  accomplished  by  forming  a 
polynomial  series  approximation  to  the  frequency  response.  By  approximating 
the  frequency  response  with  a  finite  series  of  terms#  the  h^  can  be  computed 
by  direct  calculation.  The  Remez  exchange  method^»14  perhaps  the  most 
commonly  used  algorithm  for  generating  a  finite  series  approximation  of  the 
frequency  response  function. 

As  discussed  previously#  the  passive  FTIR  application  requires  the 
development  of  narrow-bandpass  filters  that  can  Isolate  specific  spectral 
bands.  Unfortunately*  as  the  bandpass  narrows#  It  takes  Increasingly  more 
terms  to  approximate  the  frequency  response.  This  causes  the  corresponding 
time-domain  filter  to  become  Increasing  large  In  terms  of  the  number  of  h^ 
required.  Filters  generated  for  the  passive  FTIR  application  by  this 
conventional  approach  often  required  as  many  as  100  terms.  Even  If  a  100-term 
filter  Is  only  applied  to  100  Interferogram  points#  10,000  multiplications  are 
required.  This  results  In  no  savings  over  applying  the  FFT  to  the  entire 
Interferogram.  For  this  reason*  it  was  Judged  essential  that  methodology  be 
developed  for  reducing  the  number  of  required  filter  terms. 

Regress Jan  Approach  ta  £IR  Elite, r  Generation.  As  noted*  eq.  3  is  a 
linear  model  that  relates  a  filtered  Interferogram  point#  y^#  to  a  series  of 
raw  Interferogram  points#  x^._n.  Regression  analysis  Is  the  branch  of 
statistics  that  focuses  on  techniques  for  the  construction  and  evaluation  of 
such  models.  In  regression  terms#  the  y^  define  a  dependent  variable#  while 
the  x^-»n  define  a  set  of  Independent  variables.  The  h^  are  regression 
coefficients. 

It  was  Judged  potentially  fruitful  to  Investigate  the  application  of 
regression  techniques  to  the  construction  of  FIR  filters.  Two  motivating 
factors  led  to  this  conclusion.  First#  In  any  regression  model#  some 
Independent  variables  are  always  more  statistically  significant  than  others. 
These  significant  variables  contribute  most  to  the  explanation  of  variance  In 
the  dependent  variable.  It  was  hypothesized  that  If  regression  procedures 
could  be  used  to  estimate  the  h|<#  It  may  be  possible  to  delete  some  of  the  x^_n 
terms  as  being  Insignificant.  The  resulting  filters  would  thus  have  fewer 
terms.  Second#  an  Inspection  of  eq.  I  reveals  that#  In  a  conventional  FIR 
filter#  the  h^  depend  only  on  the  frequency  response  function.  In  one  sense# 
this  is  an  advantage  In  that  a  set  of  h^  can  be  used  to  filter  any  type  of 
time-domain  signal  and  any  segment  of  that  signal.  For  the  present 
application#  however#  the  only  time-domain  signal  encountered  Is  the 
Interferogram  produced  by  the  Interferometer  and  detector  of  the  passive  FTIR 
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spectrometer.  Moreover*  as  demonstrated  previously*  It  should  be  possible  to 
filter  only  a  short  segment  of  that  Interferogram.  It  was  hypothesized  that 
fewer  filter  terms  would  be  required  If  a  regression  approach  could  be  used  to 
derive  a  set  of  h|<  that  are  dedicated  to  filtering  specific  segments  of 
passive  FTIR  Interferograms. 

Generation  of  FIR  Filters  by  Stepwise  Regression.  The  generation  of  an 
FIR  filter  via  regression  methodology  begins  as  before  with  the  definition  of 
the  frequency  response  of  the  filter.  For  example*  to  generate  an  SFg  filter* 
a  Gaussian  function  would  typically  be  defined*  centered  at  the  SFg  peak 
maximum  and  possessing  a  certain  width.  The  exact  position  and  width  of  the 
bandpass  define  parameters  to  be  optimized.  As  noted  previously*  If  this 
Gaussian  function  Is  multiplied  by  a  sample  spectrum*  a  filtered  spectrum 
results.  The  Inverse  Fourier  transform  of  this  filtered  spectrum  Is  a 
corresponding  filtered  Interferogram.  If  the  desired  time-domain  filter  were 
available*  the  action  of  the  filter  on  the  sample  Interferogram  would  produce 
a  filtered  Interferogram  Identical  to  that  obtained  through  the  Fourier 
transform/Gaussian  multlpl Icatlon/ Inverse  transform  step  outlined  above. 

Thus*  points  In  this  generated  filtered  Interferogram  define  the  dependent 
variable  for  the  regression  analysis.  The  specific  points  used  would  be  those 
defining  the  Interferogram  segment  for  which  the  dedicated  filter  Is  desired. 
The  selection  of  which  points  to  use  reflects  an  additional  parameter  to  be 
optimized. 

A  standard  multiple  linear  regression  analysis  can  thus  be  performed  to 
derive  the  desired  h^.  In  this  calculation*  a  pool  of  potential  Independent 
variables  can  be  used  at  k=0,  ...*  nmax.  These  variables  would  be  derived 
from  points  In  the  same  raw  Interferogram  whose  transformed  spectrum  was  used 
In  the  generation  of  the  dependent  variable.  The  specific  Interferogram 
selected  for  the  computation  does  not  appear  to  matter*  although  It  has  been 
common  practice  throughout  this  work  to  base  the  calculation  on  an 
Interferogram  whose  spectrum  contains  the  spectral  band  that  Is  serving  as  the 
target  of  the  filter  development. 

A  stepwise  regression  algorithm^  is  used  to  select  those  x^-k  that 
contribute  most  to  the  model.  This  algorithm  begins  by  selecting  the  single 
Independent  variable  that  has  the  highest  correlation  with  the  dependent 
variable.  A  one-term  model  Is  then  formed  based  on  this  selected  variable. 
Consecutive  terms  are  added  to  the  model  In  a  stepwise  manner.  At  each  step* 
the  variables  remaining  In  the  pool  are  evaluated  for  their  correlation  with 
the  variance  In  the  dependent  variable  that  has  not  been  explained  by  the 
terms  previously  selected.  The  variable  chosen  though  this  procedure  Is  added 
to  the  model*  and  the  process  Is  repeated  until  no  remaining  variables  In  the 
pool  meet  a  minimum  standard  of  correlation.  This  correlation  test  Is 
typically  referenced  to  a  statistical  £-d1str1but1on. 

As  an  example*  an  SFg  filter  for  Interferogram  points  175-250  was 
generated  based  on  a  Gaussian  frequency  response  centered  at  941.1  cm"*  and 
possessing  a  FWhW  of  54,0  cm-*.  An  Interferogram  collected  during  a  ground- 
based  data  run  was  used  In  the  regression  procedure.  The  Interferogram  was 
normalized  based  on  the  sum  of  squares  of  points  175-250.  This  normalization 
method*  based  on  the  selected  segment*  Is  the  best  method  found  for  correcting 
for  differences  In  overall  Infrared  energy  between  Interferograms.  For  the 
stepwise  regression*  a  pool  of  variables  was  used  from  k»0  to  70.  To  be 
selected*  variables  had  to  be  significant  at  the  95%  level,  based  on  the  £- 
distribution.  The  regression  procedure  yielded  a  20-term  filter.  The 
correlation  coefficient  for  the  regression  was  0.992*  corresponding  to  98.39% 
(100  times  0.992  squared)  of  the  variance  In  the  dependent  variable  being 
explained  by  the  derived  model. 

Figure  12  presents  the  results  of  applying  this  filter  to  the  six  raw 
Interferogram  segments  depicted  previously  In  Figure  5.  Interferogram 
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Figure  12.  Points  175-250  In  the  six  Interferogram  segments  of  Figure  5  after 
applying  the  20-term  filter  computed  by  stepwise  regression. 


Figure  13.  Sum  of  squares  after  applying  the  20-term  filter  vs.  Interferogram 
number  for  each  Interferogram  In  a  typical  helicopter  data  run. 
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Figure  14.  Sum  of  squares  after  applying  the  20-term  filter  vs.  Interferogram 
number  for  each  Interferogram  In  a  typical  ground  vehicle  data  run. 

segments  3,  5,  and  6,  thought  to  contain  SFg  Information,  are  now  clearly 
different  from  the  remaining  segments,  thought  to  contain  no  SFg  Information. 

Figures  13  and  14  result  from  applying  this  filter  to  each  Interferogram 
In  two  experimental  runs.  The  Interferogram  used  In  the  generation  of  the 
filter  was  not  taken  from  either  of  these  two  runs.  Before  applying  the 
filter,  each  Interferogram  was  normalized  based  on  the  sum  of  squares  of 
points  175-250.  After  filtering  points  175-250  In  each  Interferogram,  the  sum 
of  squares  was  computed  across  the  range  of  filtered  points.  In  the  figures, 
the  computed  sums  of  squares  are  plotted  vs.  Interferogram  number.  As 
motivated  by  the  discussion  of  Figures  9-11,  the  peaks  in  Figures  13  and  14 
correspond  to  Interferograms  In  which  SFg  Information  Is  present.  The 
Interferograms  used  In  the  generation  of  Figure  13  were  collected  with  the 
spectrometer  mounted  on  a  helicopter  flying  at  80  knots  and  at  an  altitude  of 
1000  feet.  The  spectrometer  was  mounted  on  a  ground  vehicle  moving  at  10  mph 
during  collection  of  the  data  depicted  in  Figure  14.  Both  data  sets  used 
represent  worse-than-average  helicopter  and  ground-based  runs  In  terms  of  the 
strength  of  the  SFg  signals  present.  Much  higher  slgnal-to-nolse  Is  always 
observed  In  the  helicopter  runs,  as  the  spectrometer  Is  looking  at  the  target 
gas  against  the  earth,  a  strong  Infrared  radiator. 

The  center  of  the  filter  bandpass,  the  filter  width,  and  the 
Interferogram  segment  used  are  parameters  that  must  be  optimized  for  each 
spectral  band  Investigated.  The  parameter  values  used  above  were  selected 
based  on  evaluations  of  what  we  term  ideal  filter  responses.  Across  a  data 
run,  this  procedure  involves  Fourier  transforming  each  Interferogram, 
multiplying  the  computed  spectrum  by  the  desired  filter  frequency  response, 
and  then  Inverse  transforming  the  filtered  spectrum  back  to  the  Interferogram 
domain.  The  sum  of  squares  Is  then  computed  on  this  Ideally  filtered 
Interferogram.  A  plot  of  the  sums  of  squares  vs.  Interferogram  number  can 
then  be  evaluated  In  terms  of  slgnal-to-nolse.  This  effectively  evaluates  the 
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degree  of  optimization  of  the  filter  frequency  response  and  the  selected 
Interferogram  segment.  The  frequency  response  and  Interferogram  segment  are 
sought  that  maximize  this  slgnal-to-nolse  ratio. 

This  procedure  assumes*  of  course*  that  a  time-domain  version  of  the 
filter  can  be  generated.  Figure  15  depicts  the  Ideal  filter  response  across 
the  helicopter  data  set  described  above.  Interferogram  normalization  based  on 
the  sum  of  squares  of  points  175-250  was  used  as  before.  The  frequency 
response  and  Interferogram  segment  used  were  Identical  to  that  employed  In  the 
generation  of  the  20-term  filter.  It  Is  clearly  evident  from  an  Inspection  of 
the  figure  that  the  computed  time-domain  filter  djes  not  duplicate  the 
performance  of  Its  frequency-domain  counterpart. 

The  limitations  of  the  stepwise  regression  procedure  were  further 
Illustrated  when  attempts  were  made  to  construct  a  multiple-bandpass  filter 
for  DMMP.  Very  low  correlation  coefficients  were  observed  In  the  regression 
analysis. 

FIR  Matrix  Filters.  An  Inspection  of  eq.  1  reveals  that  each  point  on 
the  filtered  Interferogram  Is  derived  through  the  separate  evaluation  of  the 
convolution  Integral  for  that  point.  As  noted  previously*  a  time-domain 
filter  must  approximate  this  convolution  Integral  at  each  point.  In  a 
conventional  FIR  filter*  however*  the  same  filter  coefficients  and  the  same 
surrounding  points  are  used  to  define  each  filtered  point.  To  achieve  the 
best  approximation  of  the  convolution  Integral*  however*  It  can  be  argued  that 
a  different  filter  model  should  be  used  at  each  point.  For  each  point*  such  a 
filter  would  have  a  different  set  of  filter  coefficients  and  be  based  on  a 
different  set  of  surrounding  points. 


Figure  15.  Sum  of  squares  vs.  Interferogram  number  for  the  Ideal  response  of 
a  Gaussian  filter  centered  at  941.1  cm"*  and  having  a  FWHM  of  54.0  cm"l. 
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By  deriving  separate  filters  to  be  used  on  each  Interferogram  point#  It 
was  hypothesized  that  the  overall  accuracy  of  the  filtering  process  could  be 
Increased.  To  perform  the  derivation  of  such  a  filter#  however#  a  set  of 
Interferograms  must  be  employed.  For  example#  a  filter  to  be  used  on 
Interferogram  point  175  must  be  derived  through  the  use  of  point  175  across  a 
set  of  Interferograms.  It  Is  clear  that  the  use  of  a  set  of  filters  may  offer 
Increased  accuracy.  It  Is  also  true  that  such  a  procedure  may  be  faster 
computationally.  Since  each  filter  Is  now  performing  a  simpler  task#  fewer 
terms  In  each  filter  may  be  needed. 

The  derivation  of  a  set  of  filters  Is  a  direct  extension  of  the  stepwise 
regression  methodology  discussed  previously.  Separate  stepwise  regressions 
are  performed#  one  for  each  Interferogram  point  treated.  The  result  of  this 
process  Is  a  set  of  filters  stored  as  a  matrix  of  filter  coefficients.  For 
this  reason#  we  have  termed  the  resultant  filter  set  a  "matrix  filter." 

This  methodology  was  tested  by  developing  a  matrix  filter  for  SF5  based 
on  the  same  frequency  response  and  Interferogram  segment  used  above.  Two 
hundred  Interferograms  from  a  ground-based  data  run  were  used  In  this 
calculation.  Interferogram  normalization  based  on  points  175-250  was  used  as 
before.  The  pool  of  Independent  variables  used  at  each  point  was  defined  as 
k=0  to  70.  To  be  judged  significant#  Independent  variables  had  to  meet  a  97% 
£-cr1ter1on.  The  97%  criterion  was  used  after  discovering  that  a  95%  £ 
allowed  more  terms  Into  the  models  than  were  necessary.  No  performance  drop¬ 
off  was  observed  In  switching  to  a  97%  cutoff.  Subsequent  re- Investigation  of 
the  20-term  standard  FIR  filter  revealed  that  a  95%  £  cutoff  was  necessary  In 
that  case  to  maintain  filter  performance. 

In  deriving  the  matrix  filter  coefficients#  It  was  discovered  that  when 
computations  are  carried  out  across  a  set  of  Interferograms,  the  presence  of 
random  low-frequency  noise  prevents  the  regression  procedure  from  working 
properly.  For  this  reason#  each  Interferogram  was  preprocessed  with  a  4-term 
low-frequency  cutoff  filter  to  eliminate  low-frequency  noise. 

Figures  16  and  17  depict  the  results  of  the  regression  procedure.  Figure 
16  Is  a  plot  of  the  percentage  of  variance  In  the  dependent  variable  explained 
by  each  filter  model  vs.  Interferogram  point  number,  while  Figure  17  Is  a  plot 
of  the  number  of  filter  terms  required  vs.  Interferogram  point.  With  only  a 
few  exceptions,  the  percentage  of  variance  explained  Is  high.  Fourteen  filter 
terms  are  needed  on  average#  with  the  largest  filter  requiring  21  terms.  The 
plot  In  Figure  16  reveals  a  periodic  phenomenon  that  Is  presumably  related  to 
the  filter  bandpass.  As  depicted  In  Figure  12#  those  points  on  the  waveform 
passed  by  the  filter  very  nearly  define  a  single  cosine  function.  The  points 
near  zero  on  this  waveform  would  be  expected  to  have  low  slgnal-to-nolse.  It 
seems  plausible  that  the  regression  computation  for  these  points  would  be 
dominated  by  random  (l.e.  uncorrelated)  noise#  thus  resulting  In  a  lower 
percentage  of  variation  explained. 

To  test  this  hypothesis#  the  computed  matrix  filter  was  applied  to  the 
three  SFg-actlve  Interferograms  discussed  previously  In  relation  to  Figures  5 
and  12.  The  percentage  of  variance  explained  In  the  filter  generation  was 
plotted  vs.  the  corresponding  Intensity  of  each  filtered  point.  This  plot  Is 
depicted  In  Figure  18.  The  general  trend  In  the  plot  Is  for  those  points  with 
low  percentages  of  variance  explained  to  be  associated  with  Intensities  near 
zero.  Correspondingly#  there  are  no  points  near  the  extrema  (0.04#  -0.04) 
that  exhibit  low  correlation  In  the  filter  generation. 

The  generated  matrix  filter  was  tested  by  applying  It  to  a  number  of 
Interferograms  not  Included  In  the  regression  computation.  As  an  Initial 
test#  the  filter  was  applied  to  the  same  six  Interferograms  described 
previously  In  the  discussion  of  Figures  5  and  12.  Figure  19  Is  a  plot  of 
these  Interferogram  after  filtering.  The  SF5-act1ve  Interferograms  stand  out 
somewhat  more  cleanly  In  this  plot  than  In  Figure  12. 
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Figure  17.  Plot  of  number  of  filter  coefficients  required  vs.  Interferogram 
point. 
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Figure  19,  Plot  of  points  175-250  In  six  Interferogram  segments  after 
application  of  the  matrix  filter.  Counting  from  the  top,  the  third,  fifth, 
and  sixth  Interferogram  segments  are  thought  to  contain  SFg  Information. 
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The  filter  was  next  applied  to  the  helicopter  and  ground-based  data  runs 
discussed  In  relation  to  Figures  13  and  14.  The  Interferograms  used  In  the 
generation  of  the  filter  were  not  part  of  either  of  these  data  sets.  The  4- 
term  low-frequency  cutoff  filter  was  applied  to  each  Interferogram#  and  the 
normalization  procedure  based  on  points  175-250  was  applied.  Following 
application  of  the  filter#  the  sum  of  squares  was  computed  across  points  175- 
250.  Figure  20  is  a  plot  of  the  computed  sum  of  squares  vs.  Interferogram 
number  for  the  helicopter  run#  and  Figure  21  Is  the  corresponding  plot  for  the 
ground-based  run.  Figures  20  and  21  can  be  compared  directly  to  Figures  13 
and  14,  respectively.  From  this  comparison#  It  Is  clear  that  the  matrix 
filter  produces  superior  results  In  terms  of  slgnal-to-nolse. 

An  additional  advantage  of  the  matrix  filter  approach  Is  found  In  the 
polnt-by-polnt  feedback  from  the  regression  procedure.  Figure  16  Indicates 
that  some  points  can  be  filtered  more  reliably  than  others#  based  on  the 
success  of  the  filter  generation  procedure.  An  experiment  was  conducted  In 
which  the  25  points  below  the  90*  level  in  Figure  16  were  eliminated  from  the 
sum  of  squares  calculation.  Figures  22  and  23  are  the  resulting  sum  of 
squares  plots  based  on  the  remaining  51  points  In  the  175-250  segment.  Figure 
22  depicts  the  helicopter  run#  while  Figure  23  depicts  the  ground-based  run. 
The  slgnal-to-nolse  for  these  plots  appears  slightly  better  than  that  observed 
for  Figures  20  and  21#  validating  the  point  deletion  procedure. 

The  matrix  filter  approach  was  also  evaluated  for  Its  ability  to  generate 
multiple  bandpass  filters.  As  noted  previously#  DMMP  has  three  characteristic 
spectral  bands.  Based  on  the  Ideal  filter  response  procedure  outlined  for 
SFq,  the  three-band  frequency  response  depicted  In  Figure  24  was  Judged 
optimum.  The  Ideal  response  for  this  filter  Is  depicted  In  Figure  25.  Based 
on  an  Inspection  of  spectra#  Interferograms  277-360  are  thought  to  contain 
DMMP  signals.  Interferogram  points  175-250  were  again  found  to  be  optimum. 

One  modification  to  the  previous  procedure  was  required#  however.  In  the 
Investigation  of  the  DMMP  Ideal  response#  Interferogram  normalization  based  on 
the  entire  interferogram  (points  1-1024)  was  found  to  produce  markedly  better 
results  than  normalization  based  solely  on  points  175-250.  The  reason  for 
this  discrepancy  with  the  SFg  results  Is  unclear.  Unfortunately,  only  a  very 
limited  amount  of  passive  FTIR  DMMP  data  Is  available  to  help  pinpoint  this 
phenomenon.  This  problem  Is  made  worse  by  very  low  slgnal-to-nolse  ratios  In 
the  available  DMP  Interferograms. 

The  matrix  filter  for  DWP  was  computed  based  on  200  Interferograms.  The 
Interferograms  were  normalized  based  on  points  1-1024#  and  the  4-term  low- 
frequency  cutoff  filter  was  applied  as  before.  The  pool  of  Independent 
variables  used  was  defined  by  k=0  to  70.  Variables  had  to  meet  a  95*  £ 
criterion  In  order  to  enter  the  filter  models.  The  computed  filter  exhibited 
significantly  lower  values  for  the  percentage  of  variance  explained  In  the 
dependent  variable  than  the  filter  computed  previously  for  SFg.  Figure  26  Is 
a  plot  of  percent  variance  explained  vs.  Interferogram  point#  while  Figure  27 
plots  the  number  of  filter  terms  required  vs.  Interferogram  point.  A  more 
complex  structure  Is  seen  In  Figure  26  than  was  seen  In  the  corresponding  plot 
for  the  SFg  filter#  although  a  definite  pattern  exists.  The  filters  average 
nine  terms.  These  models  are  smaller  than  those  computed  for  SFg#  as  fewer 
terms  In  the  DMMP  calculation  met  the  £-cr1ter1on.  Relaxing  the  £-cr1ter1on 
for  the  entry  of  variables  Increases  the  sizes  of  the  models#  but  the 
percentage  of  variance  explained  does  not  Increase  significantly  and  the 
larger  models  do  not  exhibit  better  performance.  We  feel  the  low  slgnal-to- 
nolse  of  the  DMMP  Interferograms  Is  being  manifested  In  the  overall  lower 
correlations  and  the  smaller  number  of  statistically  significant  variables. 

Following  the  procedure  found  useful  for  the  SFg  matrix  filter#  those 
points  exhibiting  low  percentages  In  Figure  26  were  not  Included  In  the  sum  of 
squares  calculation  after  application  of  the  filter.  An  80*  cutoff  was  chosen 
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Figure  20 
filter  vs 


Figure  21. 
filter  vs. 


Sum  of  squares  over  points  175-250  after  applying  the  matrix 
Interferogram  number  for  the  helicopter-based  data  run. 


Sum  of  squares  over  points  175-250  after  applying  the  matrix 
Interferogram  number  for  the  ground-based  data  run. 
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Figure  22.  Sum  of  squares  over  51  points  exceeding  the  90%  variance  cutoff 
vs.  Interferogram  number  for  the  helicopter-based  data  run. 


Figure  23.  Sum  of  squares  over  51  points  exceeding  the  90%  variance  cutoff 
vs.  Interferogram  number  for  the  ground-based  data  run. 
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Figure  26.  Plot  of  percentage  of  variance  In  the  dependent  variable  explained 
vs.  Interferogram  point  for  the  DMKP  matrix  filter. 
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Figure  27.  Plot  of  number  of  filter  terms  vs.  Interferogram  point  for  the 
OMMP  matrix  filter. 


Figure  28.  Sum  of  squares  over  47  points  exceeding  the  80JE  variance  cutoff 
vs.  Interferogram  number  for  the  DMMP  data  run. 

arbitrarily  after  an  Inspection  of  the  figure.  Forty-seven  points  remained 
after  applying  this  cutoff.  Figure  28  Is  a  plot  of  the  sum  of  squares  over 
these  points  after  filtering  vs.  Interferogram  number  for  a  ground-based  DMMP 
data  run.  Due  to  the  small  quantity  of  DMW*  data  available,  the 
1 n ter ferog rams  used  In  the  filter  generation  were  derived  from  this  data  set. 
The  overall  slgnal-to-nolse  ratio  of  this  plot  Is  greater  than  that  of  the 
Ideal  plot.  We  believe  this  Is  due  tc  the  Identification  and  deletion  of  low 
slgnal-to-nolse  points  from  the  sum  of  squares  calculation. 

Quantttallvfl  Comparison  af  .SEg  Filter  Performance.  The  filtering  results 
presented  above  are  best  compared  In  a  quantitative  manner.  This  comparison 
was  made  based  on  two  criteria:  (1)  slgnal-to-nolse  ratio  of  the  sum  of 
squares  plots  and  (2)  computational  efficiency.  Based  on  Inspections  of 
transformed  spectra,  SFg-actlve  regions  were  established  in  both  the 
helicopter  and  ground-based  data  runs.  A  baseline  set  of  200  Interferograms 
was  selected  In  each  data  set  from  among  those  Interferograms  judged  to 
contain  no  SFg.  For  each  plot,  the  sum  of  squares  values  for  these  baseline 
interferograms  were  used  to  compute  a  standard  RMS  noise  level.  In  addition, 
simple  linear  regression  was  used  to  define  a  baseline  model  for  each  plot. 

For  each  Interferogram  In  the  SFg-actlve  regions,  the  baseline  model  was  used 
to  compute  a  predicted  baseline  value,  and  the  difference  was  computed  between 
the  sum  of  squares  value  for  the  interferogram  and  the  baseline  value.  The 
ratio  of  this  difference  to  the  RMS  noise  was  taken  as  the  slgnal-to-nolse 
ratio.  This  computation  corrects  for  sloping  baselines. 

Table  I  is  a  slgnal-to-nolse  comparison  of  the  plots  produced  by  the 
standard  and  matrix  FIR  filters  for  the  SFg  helicopter-based  data  run.  The 
slgnal-to-nolse  computation  described  above  was  used  In  the  generation  of  the 
table  entries.  Five  SFg-actlve  regions  were  found  In  this  data  set, 
corresponding  to  Interferograms  231-246,  320-333,  419-427,  532-535,  and  599- 
605.  Table  II  presents  the  corresponding  Information  for  the  SFg  ground-based 
run.  One  SFg-actlve  region  was  found,  corresponding  to  Interferograms  358- 
387. 
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Table  I 


Slgnal-to-Nolse  Comparison  of  Filter  Performance 
SFg  Helicopter-Based  Data 


(Based  on  Sum  of  Squares  After  Filtering) 


Filter 

Signal 

231-246 

-to-Nolse  Ratio  for  Interferogram  Region 

320-333  419-427  532-535  599-605 

Ideal  response 

73.48 

79.30 

57.39 

47.62 

68.02 

20-term  standard  FIR 

34.06 

36.33 

26.46 

23.68 

30.37 

Matrix  FIR,  76  points 

61.41 

59.25 

52.66 

38.85 

69.32 

Matrix  FIR,  51  points 
based  on  90%  cutoff 

81.68 

78.47 

69.53 

49.39 

90.82 

Ideal  response 

11.34 

20-term  standard  FIR 

4.69 

Matrix  FIR,  76  points 

19.55 

Matrix  FIR,  51  points 
based  on  90%  cutoff 

22.90 
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An  Inspection  of  Tables  I  and  II  serves  to  reinforce  the  conclusions 
reached  from  a  visual  Inspection  of  the  filtering  results.  The  matrix  FIR 
filter  significantly  outperforms  the  standard  FIR  filter.  This  Is 
particularly  evident  In  the  low  slgnal-to-nolse  data  found  In  the  ground-based 
data  run.  Furthermore*  the  deletion  of  poorly  modeled  points  from  the  sum  of 
squares  calculation  significantly  enhances  the  slgnal-to-nolse  ratios  of  the 
plots.  The  matrix  filter  methodology  Is  thus  judged  superior  In  terms  of  Its 
ability  to  extract  low-signal  spectral  information. 

One  unexplained  phenomenon  Is  seen  In  the  matrix  filter  producing  higher 
slgnal-to-nolse  than  the  corresponding  Ideal  response  for  the  ground-based 
data.  This  can  be  explained  when  points  are  deleted  from  the  matrix  filter 
sum  of  squares  calculation  based  on  the  argument  that  the  deleted  points 
contain  little  signal  Information.  For  the  ground-based  data*  however*  the 
matrix  filter  response  over  points  175-250  has  higher  slgnal-to-nolse  than  the 
Ideal  response  over  the  same  points.  This  phenomenon  remains  under 
Investigation. 

Table  III  compares  the  standard  and  matrix  filter  approaches  to  the  FFT 
In  terms  of  the  number  of  multiplications  required.  Also  Indicated  Is  the 
percentage  reduction  In  number  of  multiplications  gained  when  each  of  the 
filter  approaches  Is  used  In  place  of  the  FFT.  The  matrix  filter  based  on  51 
points  Is  the  fastest  of  the  three  filters.  Its  use  In  place  of  the  FFT 
effectively  results  In  a  90%  savings  In  the  number  of  multiplications 
required.  These  results  further  validate  the  selection  of  the  matrix  filter 
based  on  selected  Interferogram  points  as  the  best  available  filtering 
methodology  for  the  passive  FTIR  application. 

Quantitative  Nature  of  Filter  Responses.  The  discussion  of  Figures  8-11 
Indicated  that  quantitative  Information  may  be  present  In  the  Intensity  of  the 
filtered  Interferogram  segments.  Conventional  thinking  In  passive  FTIR  has 


Table  III 

Comparison  of  Computational  Efficiency 


Operation 

Number  of 

Multlpl Icatlons 

%  Savings  Over  FFT 

1024-point  Cooley-Tukey  FFT 

10240 

0.0 

20-term  standard  FIR  on 
points  175-250 

1520 

85.2 

*SFg  Matrix  FIR  on 
points  175-250 

1394 

86.4 

*SFg  Matrix  FIR  on  51  points 
based  on  90%  cutoff 

1065 

89.6 

^Includes  multiplications  for  4-term  low-frequency  cutoff  filter. 
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held  that  quantitative  determinations  are  not  possible.  For  this  reason* 
emphasis  has  been  placed  on  the  development  of  active  Infrared  sensors  based 
on  CO2  lasers.  The  ranging  capability  of  the  laser  systems  allows  a  remote 
concent rat Ion- path  length  (CL)  determination  to  be  made  for  a  target  analyte. 

The  placement  of  a  passive  FTIR  sensor  In  an  airborne  vehicle  makes 
quantitative  determinations  feasible*  however.  If  the  Input  optics  of  the 
spectrometer  are  directed  at  the  ground*  effectively  all  of  the  Infrared 
energy  entering  the  Interferometer  arises  from  ground  emission.  For  an 
aircraft  In  level  flight*  this  translates  to  a  fixed  path  length  for  the 
spectroscopic  determination*  directly  analogous  to  a  fixed-source  laboratory 
experiment.  Moreover*  the  altimeter  of  the  aircraft  provides  a  direct  measure 
of  the  path  length. 

It  would  thus  seems  feasible  to  develop  quantitative  calibrations  for 
target  analytes.  Toward  that  end*  we  have  evaluated  the  quantitative 
Information  present  In  filtered  Interferogram  segments.  Figure  29  Is  a  plot 
of  points  175-250  In  six  Interferograms  derived  from  the  helicopter-based  data 
run  at  1000  ft.  Figure  30  depicts  the  region  of  920-970  cm"*  In  the 
transformed  spectra  arising  from  these  Interferograms.  No  SFg  can  be  detected 
In  the  topmost  spectrum*  but  Increasing  levels  of  SFg  are  present  In  the  lower 
spectra.  Figure  31  Is  a  plot  of  points  175-250  after  application  of  the  SFg 
matrix  filter.  The  Intensity  of  the  filtered  signals  Increases  from  the  top 
to  bottom  of  the  plot  In  agreement  with  the  observed  spectra. 

To  evaluate  this  apparent  quantitative  Information*  the  spectral  bands  In 
Figure  30  were  Integrated  by  estimating  the  baseline  across  each  band  and 
Integrating  between  the  baseline  and  the  band.  The  magnitude  of  each  filtered 
segment  was  computed  by  taking  the  square  root  of  the  sum  of  the  squares  of 
the  Intensities.  Figure  32  Is  a  plot  of  the  computed  magnitudes  of  the 
filtered  segments  vs.  the  computed  areas  of  the  corresponding  spectral  bands. 
The  Individual  points  are  shown*  as  well  as  the  regression  line  derived  from 
the  points.  The  plot  Is  highly  linear.  Linear  regression  reveals  that  99. OX 
of  the  variance  In  the  dependent  variable  Is  explained  by  the  derived  linear 
model.  The  computed  i-value  for  the  significance  of  the  slope  of  the  plot  is 
19.47.  A  i-value  greater  than  4.0  Is  considered  strong.  In  addition* 
significant  dynamic  range  Is  available  In  the  plot  between  the  smallest  SFg 
band  and  the  baseline  point  (no  SFg).  These  results  confirm  that  a 
quantitative  calibration  can  be  obtained  directly  from  Information  In  the 
filtered  Interferogram  segments.  When  coupled  with  an  aircraft-based 
application,  the  quantitative  determination  of  target  analytes  by  passive  FTIR 
seems  totally  feasible. 

Pattacn  Becoan.1t ton  Methodology  Based  an  Filtered  Interferogram  Segments 

Overview  gf  Methodology.  The  digital  filtering  methodology  described 
above  effectively  serves  as  a  data  preprocessing  tool  to  reduce  a  collected 
Interferogram  to  an  Interferogram  segment  that  has  been  treated  (l.e.  windowed 
and  filtered)  to  enhance  the  Information  regarding  the  presence  of  a  targeted 
spectral  band.  The  next  processing  step  for  a  passive  FTIR  sensor  Is  a 
decision-making  step:  Does  the  filtered  segment  Indicate  that  the  target 
analyte  Is  present?  If  the  analyte  Is  present*  an  alarm  Is  presumably 
triggered.  If  the  analyte  Is  not  present,  nothing  Is  done. 

Decision-making  algorithms  for  numerical  data  are  termed  pattern 
recognition  procedures.  All  such  algorithms  are  based  on  two  premises:  (1) 
that  one  or  more  specific  data  categories  or  "classes"  can  be  defined;  and  (2) 
that  data  vectors  of  representative  members  of  each  of  these  classes  are 
available.  In  general*  pattern  recognition  algorithms  use  this  known, 
representative  data  to  produce  some  type  of  mathematical  construct  (vector, 
model*  etc.)  that  yields  a  numerical  response  when  applied  to  an  Input  data 
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Figure  29.  Points  175-250  In  six  Interferograms  derived  from  the  helicopter- 
based  data  run  at  1000  ft. 


Figure  30.  The  region  from  920-970  cm”*  In  the  spectra  derived  from  the 
Interferograms  In  Figure  29.  The  SFs  band  Increases  In  area  from  the  top  to 
the  bottom  of  the  plot. 
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Figure  31.  Points  175-250  In  the  six  Interferograms  after  application  of  the 
SFg  matrix  filter. 


Figure  32.  Plot  of  the  magnitude  of  the  filtered  Interferogram  segment  vs. 
Integrated  area  of  the  SFg  spectral  band.  The  line  shown  was  derived  from  the 
data  points  by  regression  analysis. 
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vector  or  "pattern."  The  numerical  response  produced  Is  Interpreted  to 
determine  to  which  data  class  the  Input  pattern  belongs.  For  a  two-class 
pattern  recognition  problem#  this  Interpretation  effectively  yields  a  yes-no 
response.  For  a  given  target  analyte  In  the  passive  FTIR  application#  two 
data  classes  exist#  corresponding  to  target-containing  Interferograms  and 
"clean"  interferograms.  For  the  work  reported  here#  the  Input  patterns  are 
filtered  Interferogram  segments. 

All  pattern  recognition  methods  explore  the  multidimensional  "pattern 
space"  defined  by  the  Input  patterns.  The  patterns  selected  as  representative 
of  the  data  classes  (the  "training  set")  are  used  to  characterize  the  regions 
of  this  space  "occupied"  by  the  defined  data  classes.  Two  general  approaches 
can  be  taken  to  this  Investigation. 

In  one  method#  the  characteristic  region  of  a  given  class  can  be  modeled, 
producing  a  vector  space  model  that  spans  the  region  of  the  class.  Unknown 
patterns  are  fitted  to  (l.e.  projected  onto)  the  class  models  to  determine  If 
they  lie  In  the  same  region  of  space.  The  SIMCA  method  of  Wold^  Is  an 
example  of  this  approach. 

The  second  basic  pattern  recognition  approach  attempts  to  model  the 
region  between  data  classes#  thereby  establishing  separating  boundaries. 

These  separating  boundaries  are  termed  discriminants#  and  the  overall 
procedure  Is  termed  discriminant  analysis^.  The  boundary  Is  typically 
defined  as  an  (n-l)-dlmenslonal  hyperplane,  where  n  Is  the  dimensionality  of 
the  patterns.  Mathematically#  the  hyperplane  Is  defined  by  an  n-dlmenslonal 
vector#  v#  that  Is  oriented  normal  to  the  surface  of  the  hyperplane.  The 
vector#  w#  Is  often  termed  a  "weight  vector."  For  a  two-class  problem,  only 
one  such  vector  Is  needed.  Unknown  patterns  are  assigned  to  a  data  class  by 
determining  on  which  side  of  the  boundary  they  lie.  This  calculation  Is 
defined  by 


d  =  *TX  (4) 

where  wT  Is  the  transpose  of  w#  x  Is  the  Input  pattern  vector#  and  the  scalar 
value,  d#  Is  termed  the  discriminant  score.  If  w  Is  scaled  to  unit  length# 
the  sign  of  d  determines  on  which  side  of  the  boundary  x  lies.  Iterative 
numerical  optimization  procedures  are  typically  used  to  define  w#  given  the 
Input  set  of  patterns  deemed  representative  of  the  data  classes. 

For  the  passive  FTIR  problem#  the  variety  of  possible  Infrared 
backgrounds  dictates  that  a  large  set  of  patterns  be  used  If  any  serious 
attempt  Is  to  be  made  to  define  a  pattern  space  that  Is  truly  representative. 
Thus#  the  selection  of  a  specific  pattern  recognition  approach  must  be  based 
largely  on  the  feasibility  of  the  method  for  use  with  a  large  data  set. 

Methods  based  on  the  construction  of  linear  discriminants  are  more  easily 
applied  to  a  large  data  set  than  are  the  class  modeling  techniques.  For  this 
reason#  these  methods  have  been  used  In  the  feasibility  studies  reported  here. 

Construction  of  Linear  Discriminants.  Data  from  two  ground-based  and 
three  helicopter-based  experimental  runs  were  assembled  to  define  a 
representative  set  of  Interferograms  for  the  pattern  recognition  feasibility 
study.  A  total  of  449  SFg-active  and  551  Inactive  Interferograms  were 
selected#  defining  a  total  data  set  of  1000  Interferograms.  These 
Interferograms  were  categorized  as  either  SFg-active  or  Inactive  based  on 
Inspections  of  the  corresponding  transformed  spectra.  Interferograms  were 
sought  with  SFg  levels  ranging  from  very  strong  to  the  limit  of  detection. 
Since  some  Interferograms  with  very  weak  SFg  signals  were  sought#  some 
discrepancies  may  exist  In  the  actual  classification  of  the  1000 
interferograms  as  active  or  Inactive. 

For  the  feasibility  study  reported  here#  the  pattern  vectors  used  were 
defined  In  a  straightforward  manner.  Points  175-250  in  the  Interferograms 
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were  filtered  by  use  of  the  SFg  matrix  filter  described  previously.  The 
filtered  segments  could  not  be  used  directly  for  pattern  recognition,  as  the 
Information  pertaining  to  absorption  and  emission  bands  Is  180  degrees  out  of 
phase  (see  discussion  of  Figures  9-11).  This  results  In  sign  differences  for 
the  Individual  filtered  points.  For  this  reason,  two  data  transforms  for 
removing  sign  Information  were  tested.  As  a  first  test,  the  values  In  the 
filtered  segments  were  squared,  employing  a  strategy  similar  to  that  used  In 
the  sum-of-squares  vs.  Interferogram  number  plots  presented  previously. 

Second,  the  absolute  value  was  taken  across  the  filtered  segment.  In 
comparing  the  results  of  these  transforms,  the  absolute  value  transform  was 
found  to  yield  the  best  results  In  the  pattern  recognition  trials.  It  Is 
hypothesized  that  this  transform  Is  superior,  as  It  Is  a  linear  transform.  As 
Indicated  by  eq.  4,  the  discriminants  developed  here  are  also  linear.  The 
nonlinear  effects  of  the  squaring  operation  do  not  seem  to  be  compatible  with 
the  use  of  linear  discriminants. 

The  1000  76-point  transformed  pattern  vectors  were  submitted  to  a 
standard  discriminant  optimization  procedure.  This  Iterative  procedure 
produced  a  weight  vector  that  could  place  981  of  the  1000  Interferograms  In 
the  correct  class  (active  vs.  Inactive).  The  19  mlsclasslf led  patterns  (9 
active,  10  Inactive)  appear  to  define  the  limit  of  detection. 

Testing  of  Computed  Discriminant.  The  computed  linear  discriminant  was 
tested  by  applying  It  to  Interferograms  In  several  experimental  runs.  In  each 
case,  a  plot  was  made  of  discriminant  score  vs.  Interferogram  number.  The 
discriminant  score  used  for  plotting  was  taken  as  the  mean  of  the  current 
computed  score  and  the  computed  score  of  the  previous  Interferogram.  This 
effectively  applies  a  two-point  moving- average  filter  to  the  discriminant 
scores.  This  procedure  decreases  the  Incidence  of  false  alarms  by 
downweighting  the  effects  of  single  interferograms  on  the  alarm  status. 

The  discriminant  was  first  applied  to  two  of  the  experimental  runs  from 
which  Interferograms  were  extracted  for  the  discriminant  development.  Figure 
33  Is  a  plot  of  the  computed  two-point  mean  discriminant  score  vs. 
Interferogram  number  for  a  ground-based  run,  while  Figure  34  Is  a 
corresponding  plot  for  a  helicopter-based  run.  Those  discriminant  scores 
exceeding  0.0  would  trigger  an  alarm  for  the  presence  of  SFg.  Each  of  the 
Interferograms  for  which  an  alarm  Is  clearly  Indicated  does  appear  to  contain 
SFg  Information. 

Figures  35  and  36,  respectively,  depict  discriminant  score  plots  for  the 
helicopter-based  and  ground-based  experimental  mns  used  previously  In  the 
evaluation  of  filter  performance.  None  of  these  Interferograms  were  Included 
In  either  the  filter  or  discriminant  development  work.  Again,  the  SFg-actlve 
Interferograms  are  clearly  distinguished  from  the  inactive  Interferograms. 
Although  one  false-alarm  Is  seen  In  the  helicopter-based  run  and  two  are  seen 
In  the  ground-based  run,  this  false-alarm  rate  represents  only  an  0.24X 
occurrence.  Based  on  these  results,  standard  pattern  recognition  procedures 
appear  completely  compatible  with  filtered  Interferogram  segments. 

Similar  pattern  recognition  trials  were  then  performed  with  pattern 
vectors  formed  from  the  reduced  set  of  51  points  found  during  the  matrix 
filter  development  work.  The  same  set  of  1000  Interferograms  was  used,  and 
the  absolute  value  transform  was  applied  as  before.  The  best  weight  vector 
obtainable  classified  only  850  of  the  1000  Interferograms  correctly,  however. 
Subsequent  application  of  this  weight  vector  to  the  test  data  runs  revealed 
much  poorer  discrimination  than  was  observed  with  the  76-point  discriminant. 

Results  obtained  during  the  development  of  the  matrix  filter  Indicated 
that  the  25  Interferogram  points  deleted  possessed  low  slgnal-to-nolse  for  the 
sum  of  squares  calculation.  Such  points  may  be  useful  for  discrimination, 
however.  If  the  active  or  Inactive  Interferograms  take  on  characteristic 
values  for  the  points.  This  question  was  Investigated  '  y  plotting  the 
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Figure  33.  Two-point  mean  discriminant  score  vs.  Interferogram  number  for  a 
ground-based  data  run.  Of  these  Interferograms,  392  were  used  In  the 
discriminant  development. 


Figure  34.  Two-point  mean  discriminant  score  vs.  Interferogram  number  for  a 
hel Icopter^based  data  run.  Of  these  Interferograms,  137  were  used  In  the 
discriminant  development. 
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Figure  35.  Two-point  mean  discriminant  score  vs.  Interferogram  number  for  a 
helicopter-based  data  run.  None  of  these  Interferograms  were  used  In  the 
discriminant  development. 


Figure  36.  Two-point  mean  discriminant  score  vs.  Interferogram  number  for  a 
ground-based  data  run.  None  of  these  Interferograms  were  used  In  the 
discriminant  development. 
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percentage  of  the  variance  In  the  dependent  variable  explained  In  the  matrix 
filter  computation  vs.  the  corresponding  computed  discriminant  value  for  the 
point.  This  plot  Is  depicted  In  Figure  37. 

The  discriminant  values  take  on  both  positive  and  negative  quantities* 
corresponding  to  additive  and  subtractive  terms  In  the  summation  of  the 
overall  discriminant  score  (eq.  4).  Nineteen  of  the  25  points  corresponding 
to  less  than  90%  variance  explained  have  negative  discriminant  values.  Of  the 
six  positive  discriminant  values*  five  are  less  than  0.1,  giving  them  low 
overall  weight  In  the  determination  of  the  discriminant  score.  The  points 
with  negative  discriminant  values  contribute  to  lowering  the  discriminant 
score  (l.e.  leading  to  an  Inactive  classification  of  the  Interferogram) .  In 
considering  this  point,  note  that  the  absolute  value  transform  causes  all 
values  In  the  pattern  vector  to  be  positive. 

These  results  confirm  that  for  pattern  recognition  purposes,  greater 
discrimination  between  active  and  Inactive  Interferograms  Is  found  by 
Including  all  of  the  points  In  the  Interferogram  segment.  In  this  example, 
the  Inactive  Interferograms  appear  to  have  large  values  at  the  points 
corresponding  to  low  percentage  variance  explained.  These  large  values  would 
be  expected  to  lower  slgnal-to-nolse  In  the  sum-of-squares  plot.  For  the 
pattern  recognition  procedure,  however,  they  help  to  Improve  the 
discrimination  between  active  and  Inactive  Interferograms.  This  Illustrates 
the  power  of  having  both  additive  and  subtractive  Information  In  a 
discriminant  calculation.  For  this  reason,  discriminants  typically  outperform 
simpler  procedures  such  as  applying  an  alarm  threshold  to  sum  of  squares  plots 
of  the  type  generated  during  the  filter  development  work. 


Figure  37.  Plot  of  percentage  variance  explained  In  the  matrix  filter 
calculation  vs.  computed  discriminant  value  for  Interferogram  points  175-250. 
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Critical  Evaluation  of  Time-Domain  Methodology 

The  concepts  described  above  for  time-domain  processing  of  passive  FTIR 
data  have  been  presented  to  several  groups  of  knowledgeable  FTIR 
spectroscoplsts.  Two  criticisms  of  the  approach  have  been  commonly  raised  and 
need  to  be  addressed. 

Questions  of  Spectral  Resolution.  It  has  been  argued  that  the  use  of  a 
short  Interferogram  segment  limits  the  analysis  to  only  low- resolution 
spectral  Information.  This  argument  Is  based  on  the  standard  FTIR  principle 
that  a  higher  resolution  spectrum  Is  obtained  by  collecting  a  longer 
Interferogram.  This  argument  Is  Invalid  on  two  grounds. 

First*  In  conventional  FTIR#  a  longer  Interferogram  must  be  collected  to 
obtain  high  resolution  Information  due  to  limitations  of  the  FFT  algorithm. 

The  FFT  assumes  the  Input  waveform  has  the  value  zero  outside  the  sampled 
region.  Thus,  the  algorithm  treats  a  short  Interferogram  segment  as  a 
harmonic  waveform  that  has  been  convolved  with  a  boxcar  function.  The 
transform  thus  contains  sin  x/x  components  corresponding  to  the  transform  of 
the  boxcar  function.  In  practical  terms#  this  results  In  a  seriously 
broadened  spectrum  If  the  sampled  Interferogram  segment  Is  short.  As  an 
example#  the  top  spectrum  In  Figure  38  Is  the  result  of  applying  the  FFT  to  a 
76-point  Interferogram  segment.  This  Interferogram  consisted  of  a  single 
cosine  frequency  at  941.1  cm"*.  The  transform  should  yield  a  single  spike  at 
that  frequency.  The  resultant  spectrum  Is  greatly  broadened*  however.  The 
lower  plot  In  the  figure  Is  the  spectrum  resulting  from  applying  an 
alternative  transform  technique  to  the  76-point  segment.  The  autoregressive 
transform  algorithm  of  Marple^®  attempts  to  model  the  frequencies  In  the  Input 
waveform.  No  assumption  Is  made  about  the  character  of  the  waveform  outside 
the  sampled  region.  The  resultant  spectrum  Is  very  close  to  a  single  spike. 

A  long  Interferogram  Is  required  only  when  the  FFT  Is  being  used  to  extract 
high- resolution  Information.  Arguments  based  on  Inherent  characteristics  of 
the  FFT  clearly  do  not  apply  when  the  FFT  Is  not  used. 

Second#  the  discussion  of  Figures  9-11  clearly  Indicated  that  Information 
pertaining  to  spectral  features  of  different  widths  Is  present  In  different 
concentrations  at  different  points  In  the  Interferogram.  As  shown  In  the 
figures#  the  Information  pertaining  to  the  spectral  band  of  Interest  Is 
dominant  when  other  features  (e.g.  the  overall  filter  bandpass  shape)  have 
damped  out.  The  selection  of  an  appropriate  Interferogram  segment  Is  purely 
based  on  the  spectral  signature  within  the  filter  bandpass.  Thus#  In  the  time 
domain#  questions  of  resolution  are  directed  to  location  within  the 
Interferogram#  rather  than  to  the  size  of  the  Interferogram  segment.  If  a 
single  spectral  band  Is  being  extracted  with  a  digital  filter#  the 
Interferogram  region  of  points  175-250  seems  close  to  optimum  for  the  typical 
widths  of  major  bands  In  organic  molecules. 

Effects  jjf  Phase  Errors.  Phase  errors  result  from  deviations  In  the 
sampling  points  along  the  time  axis  of  the  Interferogram.  Most  often#  these 
are  fractional  point  shifts  that  cause  the  Interferogram  to  be  slightly  out  of 
registration.  In  spectral-based  processing,  phase  errors  most  often  result  In 
small  spectral  artifacts.  A  phase  correction  algorithm  Is  typically  applied 
In  conjunction  with  the  FFT  to  remove  these  artifacts.  In  time-domain 
processing*  however*  phase  errors  are  potentially  more  damaging#  as  great 
dependence  Is  placed  on  correct  positioning  within  the  Interferogram.  The  use 
of  single-point  matrix  filters  places  further  demands  on  proper  registration 
of  the  Interferogram. 

During  the  development  work  described  In  this  report*  no  obvious 
deleterious  effects  due  to  Interferogram  phase  errors  have  been  observed.  We 
are  of  the  belief  that  most  phase  errors  encountered  are  fractional  point 
shifts  that  do  not  seem  to  have  great  effect.  An  experiment  was  conducted# 
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Figure  38.  Transformed  spectra  corresponding  to  a  76-point  segment  taken  from 
an  Interferogram  composed  of  a  single  cosine  frequency  (941.1  cm“l).  Top: 
the  spectrum  resulting  from  application  of  the  FFT.  Bottom:  the  spectrum 
resulting  from  application  of  the  autoregressive  algorithm  of  Marple^8. 

however*  to  test  the  effects  of  catastrophic  phase  errors  on  the  pattern 
recognition  results.  This  experiment  focused  on  the  helicopter-  and  ground- 
based  data  sets  whose  discriminant  score  plots  were  presented  In  Figures  35 
and  36*  respectively.  For  each  Interferogram*  a  different  point  In  the  region 
of  points  175-250  was  selected  randomly  and  deleted.  The  Interferogram  was 
then  compressed  to  simulate  a  point  being  skipped  while  sampling  the 
Interferogram.  The  matrix  filter  and  computed  discriminant  were  then  applied 
to  the  Interferogram  as  If  nothing  untoward  had  happened.  Figures  39  and  40 
are  the  resulting  discriminant  score  plots  for  the  helicopter-  and  ground- 
based  data  sets*  respectively. 

The  plots  of  discriminant  scores  are  noisier  than  before*  but  the  SFg- 
actlve  Interferograms  are  still  clearly  Indicated.  Several  additional  false 
alarms  are  now  Indicated  In  the  ground-based  plot*  but  overall*  the  serious 
phase  errors  Introduced  do  not  destroy  the  detection  capability  of  the 
algorithm.  When  It  Is  considered  that  the  majority  of  phase  errors 
encountered  are  less  severe  than  those  created  here*  It  can  be  argued  strongly 
that  the  current  time-domain  algorithm  Is  sufficiently  resistant  to  phase 
errors  to  be  usable. 

Source  of  False  Alarms.  In  a  final  evaluation  of  the  time-domain 
processing  algorithm*  Interferograms  generating  false  alarms  were  Investigated 
to  determine  the  cause  of  the  false  alarms.  Those  Interferograms  generating  a 
discriminant  score  just  above  0.0  tend  to  have  a  noisy  spectral  baseline 
across  the  filter  bandpass.  This  causes  the  filtered  Interferogram  segment  to 
be  similar  to  that  obtained  from  a  trace  level  of  SFg.  Those  false  alarms 
associated  with  discriminant  scores  significantly  greater  than  0.0  tend  to 
result  from  some  hardware  malfunction  during  collection  of  the  Interferogram. 
The  false  alarm  at  Interferogram  503  In  Figure  35  Is  one  such  example.  Figure 
41  Is  the  complete  transformed  spectrum  of  this  Interferogram.  In  processing 
this  Interferogram*  the  huge  low-frequency  noise  spike  seen  In  the  spectrum 
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Figure  39.  Two-point  mean  discriminant  score  plot  for  the  helicopter-based 
data  run  after  Introduction  of  catastrophic  phase  errors. 


Figure  40.  Two-point  mean  discriminant  score  plot  for  the  ground-based  data 
run  after  Introduction  of  catastrophic  phase  errors. 
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Figure  41.  Transformed  spectrum  of  Interferogram  503  In  Figure  35.  The  large 
low-frequency  noise  spike  leads  to  a  false  alarm. 


saturated  the  low-frequency  cutoff  filter  used  prior  to  the  matrix  filter.  A 
significant  amount  of  low  frequency  noise  was  left  In  the  Interferogram# 
ultimately  leading  to  the  false  alarm.  This  example  Illustrates  a  problem 
with  time-domain  processing.  If  the  digital  filters  In  use  break  down,  false 
alarms  may  result.  This  Interferogram  would  likely  not  give  rise  to  a  false 
alarm  In  a  spectral-based  processing  scheme.  Some  means  for  detecting  these 
catastrophic  phenomena  needs  to  be  Included  In  the  time-domain  processing  to 
prevent  false  alarms  of  this  type  from  occurring. 


CONCLUSIONS 


The  results  presented  in  this  report  confirm  the  practicality  of  an 
interferogram-based  processing  procedure  for  passive  FTIR  data.  The 
development  of  matrix  filters  represents  a  fundamentally  new  approach  to 
signal  processing  of  time-domain  data.  The  high  accuracy  of  the  filters 
allows  the  center  and  width  of  the  filter  bandpass  to  be  determined  precisely. 
In  addition#  easing  the  demands  on  Individual  filters  reduces  the  number  of 
coefficients  within  those  filters#  resulting  In  large  savings  In  computation 
time.  In  this  regard#  It  should  be  stressed  that  the  savings  In  computation 
time  reported  In  Table  III  assume  that  a  spectral-based  analysis  requires  only 
an  FFT.  In  reality#  spectral  phase  correction  would  be  used#  along  with  some 
type  of  spectral  subtraction  and/or  digital  filtering.  Therefore#  the  quoted 
savings  In  computation  time  are  conservative. 

The  significance  of  this  work  can  be  seen  In  two  areas.  First#  the 
ability  to  use  only  a  small  set  of  discrete  Interferogram  points  makes 
possible  a  new  generation  of  passive  FTIR  sensors  In  which  only  a  short 
Interferogram  segment  Is  collected.  Such  an  Instrument  could  be  built  around 
a  simpler  Interferometer  than  Is  currently  used#  thereby  Increasing  the 
reliability  of  the  sensor  and  decreasing  its  cost.  The  reduction  In 
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computational  requirements  simplifies  the  required  computer  hardware  In  an 
analogous  manner. 

The  second  point  of  significance  to  the  work  Is  the  elimination  of  the 
requirement  for  subtraction  of  a  background  spectrum.  The  two  target  analytes 
used,  SFg  and  DMMP,  were  both  detected  under  rapidly  changing  Infrared 
backgrounds.  The  filter-based  processing  strategy  should  therefore  feature 
Increased  resistance  to  background  changes  over  a  spectral  subtraction 
approach. 

In  summary,  the  research  described  here  provides  new  capabilities  for  the 
analysis  of  data  from  passive  FTIR  sensors.  When  coupled  with  recent  advances 
In  FTIR  Instrumentation,  It  should  now  be  completely  feasible  to  configure  a 
small,  economical  toxic  gas  alarm  system  for  both  stationary  and  mobile 
monitoring  appl Icatlons. 


RECOMMENDATIONS 


In  effect,  this  report  describes  a  feasibility  study  directed  at 
evaluating  the  potential  of  an  Interferogram-based  processing  scheme  for 
passive  FTIR  data.  Without  question,  the  work  Indicates  that  such  a  scheme  Is 
feasible.  For  the  algorithm  to  be  made  ready  for  actual  Implementation, 
however,  three  areas  of  further  Investigation  are  most  Important. 

First,  the  resistance  of  the  methodology  to  chemical  Interferents  must  be 
addressed.  In  this  regard,  we  refer  to  Interferents  whose  spectral  features 
lie  within  the  filter  bandpass  of  a  given  target  analyte.  Testing  of  the 
methodology  with  a  variety  of  analyte/lnterferent  combinations  Is  essential. 

Second,  a  formalized  optimization  study  must  be  performed  to  establish 
which  Interferogram  points  should  be  used  In  the  time-domain  processing.  Work 
performed  to  date  Indicates  that  the  optimum  Interferogram  segment  depends  on 
the  center  and  width  of  the  filter  bandpass.  Therefore,  for  each  analyte 
to  be  detected,  these  three  variables  must  be  optimized.  The  development  of 
formalized  procedures  for  performing  this  optimization  Is  critical  to  the 
optimization  of  the  overall  time-domain  processing  scheme. 

Lastly,  for  the  discriminant  development  work,  the  effects  of  training 
set  size  must  be  addressed.  Practical  limits  need  to  be  established  for  the 
necessary  size  of  the  training  set.  In  addition,  an  algorithm  for  the 
selection  of  training  set  members  must  be  developed.  These  considerations 
are  very  Important  If  the  discriminant  values  are  to  be  firmly  established. 

In  conclusion,  studies  such  as  those  outlined  above  will  serve  to  refine 
the  methodology.  We  are  of  firm  belief  that  the  time-domain  processing 
concepts  are  sound.  Work  must  now  proceed  to  optimize  the  specific  procedures 
used. 


47 


REFERENCES 


1.  Hanst,  P.  L.  In  "Advances  In  Environmental  Science  and  Technology"* 

Pitts*  J.  N.,  Metcalf*  R.  L.*  Eds.;  Wiley:  New  York*  1971. 

2.  Chan,  S.  H.;  Lin,  C.  C.;  Low,  M.  J.  D.  Environ.  Scl.  Techno! .  1973*  2, 

424-428. 

3.  Walter,  H .;  Flanigan*  D.  Appl .  Opt.  1975*  H»  1423-1428. 

4.  Low,  M.  J.  D. ;  Clancy*  F.  K.  Environ.  Scl.  Techno! .  1967,  1*  73-74. 

5.  Herget,  W.  F.;  Brasher,  J.  D.j  Opt.  £nfi.  1980*  12*  508-514. 

6.  Mertz*  L.  U.  S.  Dept,  of  Commerce  Rept.  No.  AD602936,  p.  120. 

7.  Block,  L.  C.;  Zachor*  A.  S.;  Appl .  Opt.  1964*  2.*  209. 

8.  Hanel,  R.  A.;  Schachman,  B.;  Clark,  F.  D.;  Prokesh,  C.  H.;  Taylor,  J.  B.; 

Watson*  W.  M. ;  Chaney*  L.  Appl .  Opt.  1970,  2*  1767. 

9.  Herget*  W.  F.  Abstracts  of  Papers*  SPIE  1981  International  Conference  on 

Fourier  Transform  Infrared  Spectroscopy,  Columbia*  SC;  pp  449-456. 

10.  Kauffman,  J.  W.  Abstracts  of  Papers*  SPIE  1981  International  Conference  on 

Fourier  Transform  Infrared  Spectroscopy,  Columbia*  SC;  pp  426-435. 

11.  Walker.  R.  A.;  Rex,  J.  D,  Ecflfi.  Sfl£.  Ehotft-QP-t.  InstEUm.  Enfl.  1979,  121* 

88. 

12.  Childers,  D.;  Durllng,  A.  Digital  Filtering  and  Signal  Processing:  West 

Publishing:  St.  Paul,  MN,  1975. 

13.  McClellan,  J.  H„;  Parks,  T.  W.  IEEE  Trans.  Circuit  Theor  1973,  QT-20, 

697-701. 

14.  McClellan,  J.  H.;  Parks,  T.  W.;  Rablner,  L.  R.  IEEE  Trans.  Audio 

Electroacoust.  1973,  AU-21,  506-525. 

15.  Draper,  N.  R.;  Smith,  H.  Applied  Regression  Analysis.  2nd  ed.;  Wlley- 

Intersclence:  New  York,  1981. 

16.  Wold,  S.;  Sjostrom,  M.  In  Chemometrlcs,  Theory  and  Practice,  B.  R. 

Kowalski,  Ed.;  American  Chemical  Society:  Washington,  DC,  1977. 

17.  Nlllson,  N.  J.  Learning  Machines:  McGraw-Hill:  New  York,  1965. 

18.  Marple,  L.  IEEE  Trans.  Acoust..  Speech.  Signal  Processing  1980,  ASSP-28. 

441-454. 


49 


